Dephasing in disordered metals with superconductive grains 
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Temperature dependence of electron dephasing time 7>(T) is calculated for a disordered metal 
with small concentration of superconductive grains. Above the macroscopic superconducting tran- 
sition line, when electrons in the metal are normal, Andreev reflection from the grains leads to a 
nearly temperature-independent contribution to the dephasing rate. In a broad temperature range 
r~ 1 (T) strongly exceeds the prediction of the classical theory of dephasing in normal disordered 
conductors, whereas magnetoresistance is dominated (in two dimensions) by the Maki-Tompson 
correction and is positive. 

PACS numbers: 74.40.+k, 72.15.Rn, 74.50.+r 



INTRODUCTION 

During last few years, a number of experimental data 
on electron transport in disordered metal films and wires 
were shown to be in disagreement with the standard 
theory of electron dephasing in normal conductors. 
Namely, at sufficiently low temperatures T < T± the de- 
phasing rate r~ 1 (T) was systematically found to deviate 
from the power-law dependence []J: 



(T/h) 3/2 T 1/2 /(k F l) 2 , 3D case, 



(T/2nhg) Hwg), 



2D case, 



(1) 



with a tendency to apparent saturation in the T — ► 
limit (g = H/e 2 R\j 3> 1 is dimensionless conductance of 
the film). Since no dephasing rate may exist at strictly 
zero temperature |2|, such a behavior indicates a pres- 
ence of some additional temperature scale(s) To (which 
may occur to be extremely low), so that in the range 
To < T < T\ the main contribution to t~ (T) comes 
from some new mechanism, different from the universal 
Nyquist noise considered in Ref. Q. Among other sug- 
gestions (the presence of localized two-level systems 0-Ej . 
nonequilibrium noise etc.) there were some specula- 
tions on a possible role of electron-electron interactions in 
t v {T) "saturation". Recent development |fj of the the- 
ory [j| have proved that perturbative account of electron- 
electron interactions does not lead to considerable cor- 
rections to Eq. (JJJ. 

In this paper we show that electron-electron interac- 
tion considered nonperturbatively can indeed be respon- 
sible for strong deviation of dephasing rate from the stan- 
dard predictions. Namely, we consider a system of small 
superconductive islands (of characteristic size a) situated 
in either bulk disordered metal matrix (3D case) or on 
the thin metal film (2D). The role of interaction here 
is to establish superconductivity in the islands, which is 
a nonperturbative effect. Such a system can exhibit |j| 
a macroscopic superconducting transition mediated fry 
the proximity Josephson coupling between the islands [S| , 



with the transition temperature T c (rii) depending on the 
concentration of the islands m. Above this transition 
electrons in the metal are normal, but Andreev reflec- 
tion of them from the superconducting islands leads to 
an additional contribution to the dephasing rate: 



l/T y (T) = lM°)(T) + l/T A (T). 



(2) 



Enhancement of dephasing rate due to superconduc- 
tive fluctuations in homogeneous systems was considered 
previously both experimentally [t| and theoretically |lfj| . 
Far above T c , the dephasing rate due to interaction in 
the Cooper channel is comparable to the dephasing rate 
h/r!p (T) due to the Coulomb interaction, being addi- 
tionally suppressed as l/ln 2 (T/T c ). Peculiarity of our 
result is that the superconductive contribution to the de- 
phasing rate in inhomogeneous systems can be the domi- 
nant one in a broad range of temperatures above T c (rii). 

To simplify calculations, we consider the model sys- 
tem |jj where superconducting (SC) islands are con- 
nected to the metal matrix via tunnel barriers with 
normal-state tunnel conductances Gt (measured in units 
of e 2 /h), which determine inter- islands resistance in the 
normal state. We are interested in the temperature range 
much below the critical temperature T c q of islands, when 
charge transport between them and the metal occurs due 
to Andreev reflection processes. We assume large An- 
dreev conductance, Ga ^> 1, thus Coulomb blocking of 
Andreev transport is suppressed. For small concentra- 
tion of the islands, rij < n c ~ exp(— itGa/^), quantum 
fluctuations destroy macroscopic superconductive coher- 
ence through the whole system even at T = 0, ITT) . 
In the opposite limit, ni 3> n c , the thermally driven 
superconductor-metal transition takes place at T c (rij) ~ 
HDn 2 J d , where D is the diffusion coefficient and d is the 
dimensionality of space. 

Here we focus on the temperature scale T ^> T c (rij), 
where macroscopic superconductivity is destroyed by 
thermal fluctuations, and the phases ipj of supercon- 
ductive order parameters on different islands fluctuate 
strongly and are uncorrelated with each other. Our main 
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FIG. 1: Schematic {m,T) phase diagram of a metal with su- 
perconducting grains. The dephasing time t£ due to Andreev 
reflection is shorter than t^, ' in a broad range above T c (m). 



result is the expression for the dephasing rate due to the 
processes of Andreev reflection from the SC islands: 



where 



G A (w) 



4. G A E C 
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3D case, 
2D case, 



(3) 



G d {uj) 



(4) 



is the (frequency-dependent) Andreev conductance of the 
island in the lowest tunneling approximation [l2T |. with 
Gp 1 = (e 2 /ft)(47TCTa) _1 for 3D spherical islands of radii a, 
and G^-(uj) = (A%g)~ l \n(D/a 2 uo), for 2D islands of radii 
a. Here a is the 3D conductivity of the metal matrix, 
g = 2hvD 3> 1 is the dimensionless film conductance per 
square, Ec — 2e 2 /C is the bare charging energy of an 
island and v is the metal density of states per one spin 
projection. 

Equation © is valid for T ^> max(T c (nj), Ec), where 
Ec oc EceT lxGA ^ is the renormalized charging energy 
(see below). In this temperature range the dephasing 
rate (|3J) is nearly temperature independent, thus exceed- 
ing the result for sufficiently small T < T*(m) ~ 
G 2 [ d (T)T c (n i ). Therefore, the window where Andreev 
reflection off the islands is the dominating dephasing 
mechanism is wide provided that G A (T) 1. 

In three dimensions we can also study the limit T <C 
Ec available at m -C n Cl where macroscopic super- 
conductivity never occurs due to quantum fluctuations. 
Here, the dephasing rate oc (T / E c ) [see Eq. l|55)l] 

vanishes at T — > in accordance with the general state- 
ment of Ref. 0- 

Below we provide brief derivation of the result J2J) and 
then discuss its physical origin and implications for ob- 
servable t v (T) in 3D and 2D systems. 



DESCRIPTION OF THE FORMALISM 

We start from the action functional S = Sd + St for 
the disordered metal (Sd) and tunnel junctions with SC 
islands (St), written in the replica form of the imaginary- 
time cr-model [3E1E3: 

S D - ^ Tr [D(VQf - At 3 EQ] , (5) 

St = -Y^S dA] Tl ' Q M Qs i ■ (6) 

Integration in Eq. © goes over the contact areas Aj , and 
7 = Gt /Aj is the tunnel conductance per unit area. The 
space- and time-dependent matrix Q(r, r, t') describing 
electron dynamics in the metal is the 4x4 matrix in the 
direct product of the spin space (subscripts a, (3 ... , and 
Pauli matrices a) and the particle-hole (PH) space (Pauli 
matrices r). In general, Q should also act in the replica 
space with the number of relicas N r — * 0. However, for 
the sake of perturbative calculations which do not involve 
closed loops of diffusive modes we can safely set N r = 1 
thus omitting the redundant replica space. The Q-matrix 
obeys the constraint Q 2 = 1 and the symmetry condition 
Q = T2Q T T2- The usual Green functions of disordered 
metal correspond to the stationary uniform saddle-point 
A of the action Sd [written in the energy representation, 
with E m = irT(2m+ 1)]: 



A a p(m, n) = 6 a p5 mn sign(E m ) r 3 . 



(7) 



Equation JSJ contains the superconductive matrix Qgj of 
the j-th island: 
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iifij (t) 



i<Pj{r) 





(8) 



Diffusion modes of the disordered metal are accounted 
for by the Q-matrix fluctuations near the saddle-point A. 
They can be parametrized as 



Q = A 



1 + W + ^W 2 + c 3 W 3 + c 4 W 4 



(9) 



in terms of the antihermitian matrix W obeying the con- 
straint {A, W} = 0, and C4 — c 3 — 1/8. In the PH space 
the matrix W is given by 



W 



d c 
- C t -d T 



(10) 



with d = —d) and c — c T describing diffuson and 
cooperon modes, respectively. These matrices acting 
in the spin and Matsubara spaces are nonzero only if 
£m£n < (diffusons) and e m e n > (cooperons). Their 
bare propagators have the form: 

2 

(d a p(m, n, r)d* a0 (m, n, r')) = — T>(r, r', s m , £„), (11a) 

nv 

2 

(c aj g(m, n, r)c* a0 (m, n, r')) = — C(r, r',e m ,e n ), (lib) 

TTV 
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where 

£>(r, r'j £m: £n) = ^(-Vn) Z>o(r, r', e m - e n ), (12a) 
C(r, r', e m , s n ) = 6»(e m e n ) X> (r, r', e m + e n ), (12b) 

and X>o(r, r', w) is the Green function of the diffusion op- 
erator: 

(-DV 2 + M) D (r, r', w) = S(r - r') (13) 

with the boundary condition V n 2?o( r 7 r', u>) = at the 
NS interface. 



In the 2D case, Ga(<jj) oc lnw which leads to an ex- 
tremely slow (In In r)correction to IIm (t) and, hence, to 
negligibly small Eq Q • To find Ec one then should take 
into account that the simple formula Q is modified in 
the lowest-frequency limit due to i) Cooper-channel re- 
pulsion in the normal metal, and ii) breakdown of the 
lowest-order tunneling approximation, both these effects 
were considered in |lll |. Below in this paper we assume 
(for the 2D case) that temperatures are not too low and 
approximation J2J is valid. 



PHASE TRANSITION 



DYNAMICS OF THE PHASE 



Thermal transition 



Integration over cooperon modes in the Gaussian ap- 
proximation yields the action functional that describes 
phase dynamics of the array [llEl- For a single island, 
this action is of the form (u)h = 2ixTk): 



4E C 



\u) k \G A {oj k ) 



(14) 

where Eq — 2e 2 /C is the bare charging energy, with 
C being the total island capacitance, and the Andreev 
conductance Ga(w) is given by Eq. (0} (here we ne- 
glect the interaction-induced corrections to G a studied 
in Refs. 0E1). 

The action H14JI had been studied extensively starting 
from the pioneering paper (cf. Ref. ^3] and references 
therein). At low enough frequencies, u> <C f^O: where 
f2o = Ga{Qo)Ec, only the second term in Eq. I|14|) is rel- 
evant and the theory becomes logarithmic provided that 
Ga{Qq) S> 1. The latter condition which prohibits the 
Coulomb blocking of tunneling will be assumed hereafter. 

Phase dynamics can be characterized by the 
imaginary-time phase autocorrelation function IIm (t) = 
l^vpir) -up(Qi)Y This correlator decays at the time scale 
h/Ec, where Ec is the renormalized effective charging 
energy, which is exponentially small in the considered 
regime of weak Coulomb blockade. For w-independent 
Ga(w) (corresponding to the 3D situation), the most de- 
tailed of existing estimates for Ec was found in Ref. 
using the two-loop renormalization group (RG) together 
with the instanton analysis: 



ttG a 



- E c (ttG a , 



(15) 



At T 3> Ec the deviation of the autocorrelation function 
Iljif(r) from 1 can be determined by means of RG; in the 
one-loop approximation [valid at IIm(t) S> 1/Ga] the 
result is Q: 



n M (r) = 1- 



4 ( G A E C 

itGa n V 2 -* 2 h T 



(16) 



The temperature T c (rii) of the thermal superconduct- 
ing transition is determined by the mean-field relation 



T c = J(T c )/2, J(T) = J2EAn,T), 



(17) 



where Ej(r, T) is the (T-dependent) energy of proximity- 
induced Josephson coupling between two SC islands at 
the distance r in d dimensions: 



Ej{r,T) = 



Grp 



ST 



(18) 



where £t — v fiD /2ttT is the thermal length, and we 
denoted Pz{x) = exp(— x) and P%{x) = K (x). Equa- 
tion (|17J) is valid if the number of relevant terms in the 
sum for J(T C ) is large, otherwise the transition is not of 
the mean-field type, but Eq. i|17fl can still serve as an 
estimate for T c . 

The nature of the transition in d dimensions is deter- 
mined by the parameter 5d- 



5,= 



3r 2 (fc F a) 4 
8vhDb 4:{k F l)(k F b) ' A{k F l){k F b)' 



G^p 



3tt 2 G 2 t 



G^p 



IvhD 4g 



(19a) 
(19b) 



which is an estimate for Ej(b, T) /T at T = hD/2-Kb 2 , and 
b = n i l ' d is the typical distance between the islands. In 
Eq. i|19a|) we expressed G F = rkpAj/iir 2 through the 
characteristic transmission coefficient r -C 1 of the S-I-N 
tunnel barrier. 

In three dimensions the parameter $3 can be arbitrary 
compared to 1. However, in two dimensions the param- 
eter 62 is bounded from below by the requirement of 
weak Coulomb blockade: G A (^lo) = (<Wtt) hx(l/di) ^> 1, 
where we estimated the island's capacity as C ~ a 2 /di, 
with di being the width of the insulating barrier. This 
condition requires ^2 ^> 1- Otherwise the transition is 
driven by quantum fluctuations and occurs at Ec ~ J ■ 



If Sd <C 1 then T c <C D/2nb 2 , the Josephson coupling 
is long-range and the mean-field equation l|17|) gives for 
the transition temperature: 



G 



16u 



In — = 



hD 



1 



, ,, Tr^ln-, <5 d « 1. (20) 



If <5d ^> 1 then T c 3> hD/2nb 2 and the Josephson cou- 
pling is short-range. The transition temperature can be 
estimated as 



T r = 



fcn 2 / d 



2tt 



Quantum transition 



S d » 1. (21) 



Quantum transition can be described within the low- 
est tunneling approximation only in three dimensions 
(cf. for discussion of quantum phase transition in 
a more complicated 2D case). Upon decreasing m, the 
transition temperature defined by Eq. H20|l lowers even- 
tually below Ec, then quantum fluctuations should be 
taken into account. At some critical concentration n c 
the temperature of the superconductive transition van- 
ishes, marking the point of a quantum phase transition. 
The point of the quantum transition is determined by 
the equation similar to (|17fl : Eq = J~(0) (cf. 7] for more 
details). However, the zero-temperature value of the in- 
tegrated Josephson proximity coupling i7(0) cannot be 
determined by the simple formula i|18|) due to logarith- 
mic divergency of the resulting expression. This diver- 
gency is cured by the account of the Cooper-channel re- 
pulsion constant in the metal A„ || leading to J{0) — 
G^rii/lQvXn. As a result, the critical concentration n c 
is found to be 



\Q-Kv\ n E t 



c 



where Eq is defined in Eq. 



COOPERON SELF-ENERGY 

In the presence of SC islands, cooperon modes are no 
longer gapless. To obtain the cooperon self-energy due 
to Andreev reflection we calculate the correction to the 
action in the lowest tunnel approximation: 



6S = 



} - (S^SP) + {bD , (23) 



where the vertices and Srp come from expansion of 
the actions JHJ) and ©, respectively, to the order W l . 
The second order in Gt a ppro ximation ill'MI) is valid pro- 
vided that Gt C Gjj |12l llq|. The corresponding dia- 
grams are shown in Fig. [21 Their sum is independent of 





FIG. 2: Diagrams for the cooperon self-energy in the second 
order over Gt- Shadowed blocks are cooperons and diffusons, 
dots denote Andreev reflections from the dot, and wavy lines 
stand for the phase correlation function Il(wfe). 



the certain form of the parametrization 0. Averaging 
in Eq. goes over phase ipj (r) dynamics and bare dif- 
fusive modes <|llll . It is important that at T ^> T c the 
phases on different islands are uncorrelated with each 
other. Upon averaging, one obtains the cooperon part 
of the induced action , which in the long- wavelength 
limit takes the form: 



SS ( 



rfr (m,n,r)\' 



(24) 



where c is the cooperon part of the matrix W, Eq. 1|1U[I . 
and 



TliG'rp 

IQv 2 



[D(r, r'; m,n — k) 



C(r,r'; 



dAdA' 
A 2 

-k)}U M (k) + {m 



-n}, (25) 



and IIm (fc) is an imaginary- frequency version of the auto- 
correlation function Hm{t). Equation l|24|) is valid pro- 
vided that the cooperon wave vector q is smaller than 
the inverse separation between the islands, q <C n\ , 
which allows to pass from the discrete sum over the is- 



lands to the uniform integration over r: ^ . 



I dr. 



The self-energy S m „ determines the \ow-q behavior of the 



(22) cooperon: C(q, m, n) — (Dq 2 



y )- 1 



Integrating diffusive modes over the area A of the con- 
tacts yields the normal-metal resistance G^ 1 which com- 
bines with G\ into the Andreev conductance. After sim- 
ple algebra we obtain for e m , e n > 0: 

^ niG A 



8;/ 



T 



,k— — rr, 



k——r 



(26) 



which is written for the case of w-independent Ga- Ana- 
lytic continuation of Eq. 126|) from e m , e n > to real fre- 
quencies, ie m — * e + iO, ie n — > e' + iO, yields the cooperon 
self-energy 



a niGA 
S(e, £ ) = — i 



$v J-oo 2vr 

+ n R {n)F(e-n) + n A (n)F(e' + n)} (27) 

in terms of the Keldysh, retarded and advanced compo- 
nents H K ' R,A (£l) of the phase correlation function, and 
F(fl) = tanh(fi/2T). 
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To study the quantum corrections to conductivity at 
zero frequency we set s' = — e leading to the cooperon 
decay rate 7(e) = — 7i -1 E(£, — e): 



lie) 



4Hv 



r°° dtt 

I -n(n){i-F(fl)P(n- E )}, (28) 



where we used the identity Xl K (tt) = -2ffl(fi), where 
II(ri) is the Fourier-transform of the real-time sym- 
metrized autocorrelation function of the island's order 
parameter H(t) — (cos[ip(t) — cp(Q)]). Another useful rep- 
resentation for 7(e) follows from Eq. I|28|) by means of 
the inverse Fourier transformation into the time domain: 



range where 1/t£ is the main contribution to the de- 
phasing rate. Using the 3D expression for G a one finds 



Tf>i)~ -x 2 J\Tk F l) 

T 



4/3 



(30) 



where Xi — (47r/3)a 3 ni is the volume fraction of the su- 
perconductive material in the matrix. From the low- 
temperature side applicability of the 3D result J2J is 
limited by the thermal transition temperature T^ D (m). 
Thus the relative width of the temperature window where 
Andreev reflection from the SC islands is the leading 
mechanism of dephasing is given by the ratio 



7(e) = — - — Tcoth — 
n ' 2hv 2T 



n(t) 



sinh^i 



dt. 



(29) 



It is interesting to note that the functional form of 
Eq. 1|29|) coincides exactly with the expression for the tun- 
neling density of states in the presence of the Coulomb 
zero-bias anomaly, cf. Eq. (58) of Ref. In the present 
case the island's phase <p(t) plays the role of the Coulomb- 
induced phase K(t) introduced in [l^, whose fluctua- 
tions give rise to the zero-bias anomaly. Then expression 
(|29|l can be rationalized with simple physical interpre- 
tation: "superconductive" contribution to the cooperon 
decay rate is just the average rate of Andreev processes 
which occur in the system. Indeed, quantum correction 
to conductivity comes from interference between differ- 
ent trajectories of the same electron; Andreev reflection 
transforms this electron into a hole, therefore destroying 
further interference. 



500 G 



3/2 



T? D (n t ) 



In 2 (n,/n ) ' 
5 0gf (np/n,) 1 / 3 
ln(noM) 



rit » n , 



rii < n , 



(31) 



where we used Eqs. i|19ajl . I|20|l and (|21|l . and defined 
ri = (8vHD/Gt) 3 such that S 3 — (n i /n ) 1 /' i . Large 
factors in Eq. I|31(l result partly from the large factor in 
Eq. (|30|l hidden there in Xi and T, and partly from writ- 
ing ln(5 3 = (1/3) ln(nj/n ). We see that the condition 
Ga 3> 1 guarantees the existence of the broad tempera- 
ture range where the dephasing time is nearly tempera- 
ture independent and given by t£. 



Lowest temperatures, T <C Ec 



DEPHASING TIME 

Now we start to analyze the consequences of the result 
B29JI . To evaluate the islands' contribution into the de- 
phasing rate, we need 7(e ~ T). Behavior of U(t ~ h/T) 
is governed by the relation between temperature T and 
the effective charging energy Ec of SC islands. 



The region of very low temperatures, T Ec, can be 
traced only at very small concentration of the island, rii < 
n c [cf. Eq. I|22|l ]. where superconductivity is absent even 
at T = due to quantum fluctuations, Here the integral 
(|29|l converges at t ~ fo/Ec and can be approximated as 



7(e) 



rijG A 
2irv 



(32) 



3D case 



Moderately high temperatures, T > Ec 



The above integral is of the order of Eq . Then the 
Andreev-reflection contribution to the dephasing rate can 
be estimated as 



At T > Ec, the integral in Eq. converges at t ~ 
h/T where H(t) is given by Eq. (|TB|) . As a result, 7(e) is 
nearly energy-independent at e ~ T and can be identified 
with the dephasing rate leading to the 3D result The 
latter is valid as long as the expression in the brackets is 
large compared to unity. 

Assuming that (T) is given by Eq. , we can es- 
timate the upper boundary T^ D {rii) of the temperature 



1 



T 



tMT) 2irhv E C 



(33) 



Since 1/r^ scales oc T it always dominates the standard 
3D result Q at very low temperatures. However, the 

crossover temperature, where = scales as nf 

and can be extremely low for small concentration of the 
islands. 
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2D case 

As explained above, staying within the lowest tunnel- 
ing approximation we can explore only the region of rel- 
atively high temperatures, (Gt / '4irg)ln(hD / 'a 2 T) -C 1, 
where fluctuations are thermal. Substituting II(i) by 1 in 
Eq. (|29|) we come to the result © for the 2D case. Here, 
contrary to the 3D case one can neglect the one-loop fluc- 
tuation correction oc lnlnT compared to the bare InT 
dependence of Ga- 

Comparing with Eq. Q one finds that the "super- 
conductive" contribution to dcphasing is dominant at 
T < T 2D ( ni ), where 



T 2D (m) =TrhDn., 



Ga(T* ) 
ln(ng) 



(34) 



The relative width of this window is then estimated by 
the ratio 



20 G A {T,) 



T 2D (n % ) ~ ln(^)ln 2 (G|/4 9 ) 



(35) 



and is large since Ga 1. 



MAGNETORESISTANCE 

Experimentally t v is determined from the magnetore- 
sistance data. For 2D systems, the low-field magnetore- 
sistance is governed by the weak localization (WL) and 
Maki-Tompson (MT) corrections which have the same 
dependence on the magnetic field [20| : 



AR{H) 
R 2 



' \a- m] Y (36) 



27T 2 ft L 



he 



with Y(x) = ]n(x)+ip(l/2 + l/x). Here a = 1 (-1/2) is 
the WL contribution in the limit of weak (strong) spin- 
orbit interaction, while the MT contribution is character- 
ized by the function (3(T) expressed through the Cooper 
channel interaction amplitude T(wk) 20] : 

= y ^(-l) m r(o; m ) - r"(^2„+i). (37) 

m ri>0 

In a uniform system far above T c , (3(T) ~ 1/ In 2 (T/T c ) 
indicating that the MT contribution is smaller but in 
general comparable to the WL contribution. 

For our system, effective attraction in the Cooper chan- 
nel emerges as a result of Andreev reflection from the 
SC islands. Formally, integration over the phases (pj(r) 
of the islands generates the standard Cooper interac- 
tion term in the action with T(ujk) = (riiG T /16u)JI(u)k)j 
where we made use of the fact that the phases of differ- 
ent islands are uncorrelated at T > T c and performed 
spacial average justified by the inequality L v 3> b. In 



the temperature range considered, II(r) is nearly con- 
stant on the time interval r € [0, 1/T], so one can use 
the static approximation H(uJk) = Sk,o/T. Substituting 
into Eq. w e obtain for T > T c 



P(T) 



i,G n 



64 vT 



(38) 



Comparing with the estimate (|34|l one finds that f3(T) ^> 
1 at T <C T*, that is magnetoresistance is mainly due to 
the MT term and thus is positive irrespectively of the 
strength of the spin-orbit scattering. 

Another relevant contribution to magnetoresistance is 
the Aslamazov-Larkin (AL) correction. In the range 
T T c , using the condition 62 3> 1, one can estimate 
A(/al < HriiD/T. Comparing with A^mt following from 
Eqs. JHEl and igEJ, one finds A^al/A^mt ~ l/<5a < 1- 
Moreover, the relevant scale of magnetic field £?al for the 
AL contribution to magnetoresistance is Bal ~ cT/eD, 
i.e., it is by factor Tt^/H 3> 1 larger than the correspond- 
ing WL scale -Bwl ~ hc/eDr v . Thus AL correction to 
magnetoresistance is much smaller than quantum (WL 
and MT) corrections, and r v can be extracted from the 
standard low-field magnetoresistance measurements. 

We believe the same conclusion to be valid in the 3D 
case. Here, however, the MT correction can be either 
large or small compared to the WL correction, depending 
on temperature and other parameters of the problem. 



DISCUSSION 

We demonstrate that small concentration of super- 
conductive islands can enhance considerably the low- 
temperature dephasing rate in disordered bulk and thin- 
film metals as seen via the low-field magnetoresistance. 
In 2D the dominant quantum correction to magnetoresis- 
tance is the Maki-Tompson one, thus magnetoresistance 
is positive in the region of interest. Throughout the whole 
range where our results are valid, Tr^/h 3> 1, which al- 
lows to neglect the energy dependence of the cooperon de- 
cay rate (|29|l . This is why magnetoresistance follows the 
standard formula l|36|) derived for uniform metal films. 

It was implicitly assumed while deriving Eq. 1|29J) that 
Lip = y Dt v is much longer than inter-island separa- 
tion b. Using Eq. J3J one finds that in 2D case for 
this condition to be fulfilled the tunnel-limit inequal- 
ity Gt/Gd < 1 is required; for 3D case the condi- 
tion L v ^> b is less restrictive. We expect that in the 
2D case with SC islands strongly (Gt >• Go) cou- 
pled to the film the "Andreev" contribution to the 
dephasing rate at moderate temperatures can be esti- 
mated analogously to Eq. ©, with the proper expres- 
sion Ga ~ Go for the Andreev conductance, leading to 
~ UiD/ ln(£r/a). Although we considered temper- 
atures much below the intrinsic transition temperature of 
SC islands T c o, our approach can be adapted for T ~ T c q. 
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We note in passing that inhomogeneous in space super- 
conductive gap function is known to affect the BSC peak 
in the density of states in a way very similar to that of 
magnetic impurities [2lj . The present results show that 
analogy between inhomogeneous superconductivity and 
magnetic impurities extends to dcphasing as well. The 
influence of the same dephasing mechanism upon other 
phase-coherent phenomena (e.g., mesoscopic fluctuations 
and persistent currents) remains to be studied. 
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